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ABSTRACT 


Micro-satellites often require the adaptation of existing propulsion systems. 
Eleetric propulsion thrusters are perhaps the best eandidates to meet these needs and ion 
engines are among the most sealable. Miniaturizing the ion engine will require novel 
eoneepts for the ionizer with perhaps novel propellants. MEMS, nanoteehnology and 
other teehnologieal advanees are expeeted to impaet on new designs. 

Our work shows that the ionization of Argon, whieh is an alternate fuel to Xenon, 
ean be aehieved at low voltages by utilizing Miero-Struetured Eleetrode (MSE) Arrays. 
Copper-elad sheets separated by a dieleetrie material (fiberglass laminate epoxy resin 
system combined with a glass fabrie substrate) of varying thiekness (0.1 mm to 0.4 mm) 
form the diseharge electrodes in the MSE arrays The wafers are drilled with an array of 
holes and this geometry serves to eoneentrate the electrie field between electrodes 
enhaneing eleetron emission at the cathode. Minimum breakdown voltages between 240 
and 280 Volts at pressures of around 100 mTorr (0.133N/m ) were eonsistently obtained 
with arrays of hole diameter ranging from 300 to 500pm. These results are eonsistent 
with eonventional Pasehen-eurves with two empirieal eonstants that arise from our 
uneonventional geometrieal arrangements and from the different material properties. 
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1. INTRODUCTION 

A. ION ENGINE HISTORY 

The theories behind ion engines have been attributed to the German scientist Dr. 
Wernher von Braun during the 1930s. At this time, the German military was more 
interested in weapons of war than interplanetary rockets, so his theories remained 
untested until the end of World War 11, when he and hundreds of fellow scientists were 
brought to the United States. Once in the US, they were able to develop their theories 
enough that in 1958, the Army Ballistic Missile Agency initiated a contract with Electro- 
Optical Systems to study ion propulsion. The result of this was a 0.1 pound-thrust engine 
developed by Hughes Research Laboratory before it stopped due to the Apollo program 
[8, 9]. 

In the early 1990s, the NASA Solar Electric Power Technology Applications 
Readiness project revived the Hughes ion thruster. This project began to study the idea of 
using a Xenon propellant within an ion engine. In 1996 one such engine was built and 
fired for over 8000 hours, making the test a success. Deep Space I was a small space craft 
that was launched onboard a small Delta 11 rocket from Cape Canaveral Air Station, 
Elorida on October 24, 1998. The propulsion system managed to propel the craft at a rate 
where it could sometimes travel at 750,000 miles per day. After the scheduled end of the 
test, the mission was extended for the opportunity to fly close to the comet Borrelly, and 
take some of the best pictures to date [8]. 

B. ION ENGINE THEORY 

Modern ion thrusters use inert gases for the propellant. The majority of these 
thrusters use Xenon, which is chemically inert, colorless, odorless, and tasteless. The 
propellant is injected upstream of the thruster and flows to the downstream end. This 
injection method is preferred because it increases the time that the propellant remains in 
the chamber thereby ensuring higher ionization efficiencies [9]. 

In a basic ion thruster, electrons are created at a hollow cathode, called the 
discharge cathode, located at the center of the engine on the upstream end. The electrons 
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flow out of the discharge cathode and are attracted to the discharge chamber walls, which 
are charged to a high positive potential by the thruster’s power supply. 

The electrons from the discharge cathode ionize the propellant by means of 
electron bombardment. High-strength magnets are placed along the discharge chamber 
walls so that as electrons approach the walls, they are redirected into the discharge 
chamber by the magnetic fields. By maximizing the length of time that electrons and 
propellant atoms remain in the discharge chamber, the chance of ionization is maximized, 
which makes the ionization process as efficient as possible at these low pressures. 

In an ion thruster, ions are accelerated by electrostatic forces. The electric fields 
used for acceleration are generated by electrodes positioned at the downstream end of the 
thruster. Each set of electrodes, called ion optics or grids, contains thousands of coaxial 
apertures. Each set of apertures acts as a lens that electrically focuses ions through the 
optics. 

Most ion thrusters use a two-electrode system, where the upstream electrode is 
charged highly positive, and the downstream electrode is charged highly negative. Since 
the ions are generated in a region of high positive charge and the accelerator grid’s 
potential is negative, the ions are attracted toward the accelerator grid and are focused out 
of the discharge chamber through the apertures, creating thousands of ion jets. The 
stream of all the ion jets together is called the ion beam. The thrust force is the change of 
momentum to the ions produced by the accelerator grid. The exhaust velocity of the ions 
in the beam is proportional to the voltage applied to the optics. 

Because the ion thruster produces a large amount of positive ions, an equal 
amount of negative charge must be expelled to keep the total charge of the exhaust beam 
neutral. A second hollow cathode called the neutralizer is located on the downstream 
perimeter of the thruster and expels the needed electrons as depicted in Eigure 1. 
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Figure 1. Ion thruster operation (From Ref. 9) 

The ion propulsion system consists of five main parts: the power source, power 
processing unit, propellant management system, the control computer, and the ion 
thruster. The ion propulsion system power source can be any source of electrical power, 
but solar and nuclear are the primary options. A solar electric propulsion system uses 
sunlight and solar cells for power generation. A nuclear electric propulsion system uses a 
nuclear heat source coupled to an electric generator. The power processing unit converts 
the electrical power generated by the power source into the power required for each 
component of the ion thruster. It generates the voltages required by the ion optics and 
discharge chamber and the high currents required for the hollow cathodes. The propellant 
management system controls the propellant flow from the propellant tank to the thruster 
and hollow cathodes [9]. 

The control computer controls and monitors system performance. The ion thruster 
then processes the propellant and power to perform work. Modem ion thmsters are 
capable of propelling a spacecraft up to 90,000 meters per second (about 200,000 miles 
per hour). The tradeoff for this high top speed is low thmst (or low acceleration) [9]. 
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Modern ion thruster units ean deliver up to 0.5 N of thrust. To eompensate for 
low thrust, the ion thruster must be operated for a long time for the spacecraft to reach its 
desired speed. Because ion thrusters use inert gas for propellant, they eliminate the risk of 
explosions associated with chemical propulsion. The usual propellant is Xenon, but other 
gases such as Krypton and Argon may be used [9]. 

C. FIELDS COVERED 

Ion propulsion is one of several propulsion methods utilized for orbital 
maintenance of spacecraft as well as interplanetary applications. With the advent of the 
Hall thruster, ion propulsion technology is beginning to take a back seat in research as 
Hall Effect technology is more readily adaptable to different applications, however, ion 
propulsion technologies should not be ignored. Current ion technology utilizes Xenon as 
the propellant of choice due to its high atomic mass, which as the mass is increased so is 
the momentum change and thus the net thrust. Unfortunately, Xenon is quite expensive. 
If more common elements such as Argon could be used the price of the propellant would 
become much cheaper [4]. 

Another problem with current ion engine technologies is that they have reached 
size limitations with current technologies. To make ion engines a viable alternative to 
other propulsion methods they must be scalable beyond where they are now. A possible 
approach to this would be the use of Micro Structured Electrode (MSE) Arrays [6]. This 
would reduce or even eliminate the ionization chamber of current engines allowing huge 
mass and size savings. It would also allow the engines to be scaled to much smaller sizes 
for applications such as micro satellite propulsion. Another benefit of utilizing MSE 
Arrays is the reduction in power requirements. Typical ion engines operate in the 
kilowatt range but by redesigning the technology to use MSE arrays the power 
requirements would drop to the range of hundreds of watts and at only a few hundred 
volts, thus requiring less power generation and a smaller spacecraft power bus. 

This thesis is concerned with utilizing MSE arrays to modify the design of the 
ionizer section. Different configurations of MSE arrays are studied in a wide pressure 
range from 10 milli-Torr to 1 Torr (or 0.122289 N/m^ to 133.289 N/m^). The 
configurations studied involved materials and design, the type of insulator and the 
structure diameters. 
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The MSB arrays were studied with respect to the electrical characteristics of the 
materials and the gas. 

The second chapter covers the test equipment setup used to for the experiments 
and offers a description of the MSB wafers. The third chapter covers ion engine theory, 
both current ion engine configurations and MSB array theory. The MSB array theory will 
cover the breakdown process at various pressures and the mechanisms for generating 
sustained discharges. The fourth chapter will cover the data and results of the 
simulations and testing of the MSB arrays. The fifth and sixth chapters will draw 
conclusions based on those results and recommendations for future studies. The thesis of 
BT Brank Perry compliments this effort [7]. 

The ultimate purpose of this research is to explore the use of MSB arrays as a 
viable alternative to the present ionizing chambers in the ion propulsion system. 
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II. EXPERIMENTAL APPARATUS 


A. VACUUM CHAMBER AND ASSOCIATED EQUIPMENT 

The primary test apparatus used in this work eonsists of a stainless 
steel/glass/plexiglass vacuum chamber, two roughing pumps, a turbomolecular vacuum 
pump, an Argon supply system, high voltage DC power supply, and various metering 
equipment such as a volt meter, ammeter and oscilloscope. The vacuum chamber is a 
cylindrically shaped glass chamber with a removable plexiglass top cover, and a bottom 
stainless steel interface where a high-vacuum gate valve connects to the turbomolecular 
vacuum pump. The stainless steel interface also has various ports to connect the two 
roughing pumps, the Argon supply system and the various electrical leads from/to the 
power supply and metering equipment. One of the roughing pumps is mounted to the 
side of the turbo pump and is used to draw the pressure of the chamber down to the tens 
of milli-Torr range before the turbo pump is started. The thesis written by Frank Perry 
covers the vacuum chamber set up in more detail [7]. For start up sequence and 
operation procedures refer to Appendix A. 

Once the Turbo molecular pump is engaged and draws the vacuum down to the 
range of 10'^ Torr the turbo pump and roughing pump are isolated and Argon is back 
filled into the chamber to the desired pressure for testing by the use of a metering valve. 
The pressure in the chamber is monitored by thermocouple sensors and a filament sensor. 

A DC power supply is used to feed voltage to the sample wafer in the chamber 
while the breakdown voltage is monitored by and oscilloscope and a voltmeter with an 
ammmeter to measure the current drawn. 

Once the breakdown voltage has been achieved the system is reset for the next 
experiment by evacuating the chamber to 10'^ torr. 

B. MICRO-STRUCTURED ELECTRODE ARRAY WAFERS 

The MSB structures used in these experiments are made from two copper layers 
sandwiching a fiberglass laminate epoxy resin system combined with a glass fabric 
substrate insulator. These Cu-dielectric-Cu structures are cut from larger sheets into 
approximately two inch by two inch pieces. These are then marked appropriately for 
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insulator thickness and nine micro-holes are drilled, using a preeision drill press, in the 
eenter in a three by three grid pattern. The holes are spaced so that they are 2 mm apart. 
Eaeh wafer is then etched using Ferric Chloride so as to strip the eopper away from the 
edges - thus leaving 10 mm of dieleetric exposed to prevent current flow across the edges 

m. 

Once etehed, they are cleaned of oils and other eontaminants by immersion into 
an aleohol bath and rubbed down. Exeess alcohol is allowed to evaporate. Each wafer is 
then inspeeted for contaminants and then plaeed into the vacuum chamber for 
experimentation as needed [7]. 
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III. THEORY OF IONIZATION ENGINE MODIFICATIONS 


A. TOWNSEND THEORY OF BREAKDOWN 

The concept behind the use of MSE arrays is to enhance the local electric field 
when a potential difference is applied to the electrodes. Within a gaseous medium such 
as Argon, the high field regions free electrons that can be amplified in an avalanche 
process. When the electric field exceeds a threshold level the number of charged 
particles multiplies exponentially and a discharge is started [6]. 

When a sufficiently strong electric field is applied, breakdown occurs. 
Breakdown is the process where a gas is converted from a non conducting material to a 
conducting one. 

From the Townsend theory of breakdown, we know that the charge carriers are 
produced by volume processes. This is depicted by the ionization coefficient a, and by 
secondary emission coefficient y. To start a self sustaining discharge, for every electron 
lost at the anode one has to be replaced by a secondary electron that was created in the 
gas or at the cathode. The ionization coefficient depicts how the electrons multiply in the 
gas medium in the direction of the electric field. The second ionization coefficient 
depicts the electron production at the cathode-gas interface. The secondary coefficient 
explains the effect of ions, photons, and neutrals [1]. 

At low pressures the production of electrons is primarily caused by ion impact on 
the cathode surface. The breakdown voltage can be shown in Equation (1) to be; 


K=- 


Bpd 


ln(pj)-l-ln( 


k) 


ln(l + y ) 


( 1 ) 


where p is the pressure of the gas medium in Torr (760 Torr = 1 atmosphere at sea level, 
1 Torr = 133.289 N/m ), d is the distance between parallel-plate anode and cathode in cm, 
A and B are constants that vary for each gas medium [1]. The breakdown voltage can be 
graphically represented by the use of a Paschen curve, which depicts the breakdown 
voltage versus pressure-times-distance as shown by plotting Equation 1 in 
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nondimensional form in Figure 2. This curve was generated using the general 
equation (Equation 2) for a Paschen curve. 


Y = 


X 

hW 


( 2 ) 


Where Y represents the Breakdown voltage and X represents the pressure times distance. 
Figure 3 shows the general Paschen curve for an Argon medium. 


Paschen Law (Y=X/ln(X)) 



On the left hand side of a Paschen curve the breakdown voltage rises rapidly as 
pressure times distance (pd) decreases due to the low possibility of ionizing collisions 
requiring a strong electric field. On the right hand side of the curve the breakdown 
voltage rises gradually as pd increases due to the probability that an electron will produce 
ionization even at lower electric fields than the left hand side. The minimum of the curve 
is where the ionization process occurs at the minimum breakdown voltage. This 
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minimum voltage region is characterized by a specific pd and these two parameters are 
the primarily concerns for the miniaturization of the ionization chamber and 
are explored in detail throughout this thesis. 
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Paschens Spark Curve for Argon 
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Figure 3. Argon Paschen Curve using a log scale. 


The limitation of the Townsend theory is its assumption of homogeneous electric 
fields and pd values. When using non-homogeneous fields and variable geometry the 
Townsend theory breaks down and it is necessary to determine the field structure that is 
being produced. Experiments with the MSE geometries differ greatly form those with 
planar geometries discussed previously. 

B. MSE ARRAY GEOMETRY 

With three-dimensional MSE array geometries it can be expected that high 

electric fields can be achieved locally within the structure. The electric field structure in 

the MSE array can be expected to be non uniform with respect to planar geometries and 
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localized within the structure. Due to the imperfect manufacturing techniques used for 
these experiments the precise geometry of the individual holes differs slightly from each 
other. The high electric field thus generated in the structure is desirable for the ionization 
of Argon producing plasma and aids in electron emission from the cathode [6]. 



Figure 4. Cross-sectional view of Cu-dielectric-Cu layers without microhole 

structures. 



Figure 5. MSB array Cross-sectional view of Cu-dielectric-Cu layers with microhole 
structures (Hole diameters from 300 pm to 500 pm). 


The three-dimensional MSB array geometry used is a three-by-three grid of holes 

in a two-inch by two-inch structure of two copper layers sandwiching a silicate insulator. 

The copper layers represent the electrodes. Figure 4 is a cross sectional representation of 
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this composite structure. Figure 5 is a eross seetional representation of the wafer with a 
hole drilled through it from top to bottom. Each of the nine holes can be a souree of a 
miero-diseharge that takes plaee between the two copper layers as they are in parallel. 
The insulation layers used for this experiment range in thiekness from 100pm to 400pm 
and the holes range in diameter from 300pm to 500pm. For further information on the 
wafer construction please refer to Perry [7]. 

C. FIELD EMISSION 

Field emission of electrons from the cathode could be a primary agent causing 
ionization in these experiments. Field emission is the proeess by whieh electrons are 
liberated from a surfaee that is under the effeet of an eleetric field. These electrons are 
then aeeelerated in the applied electrie field until they eollide with neutral speeies, in the 
ease of Argon gas these are atomic species. During the eollision, additional eleetrons are 
liberated, whieh in turn are accelerated by the electrie field and collide with other atoms 
ionizing them and produeing even more eleetrons and so on eausing an avalanche effect. 
This is an avalanche situation due to the exponential nature of the electron liberation. 

Our experiments will benefit by utilizing the field emission effeet to cause 
ionization in an Argon medium. The unique geometry of the three dimensional MSE 
arrays ereates a coneentrating effeet of an applied eleetrie field allowing the liberation of 
some eleetrons through the field emission effeet at lower energy levels than would be 
required with parallel eleetrodes. There may be ways to further enhanee their effeets 
through manufaeturing improvements and through the introduction of carbon nanotubes 
at the eathode [4]. 
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IV. EXPERIMENTAL RESULTS 


A. INTRODUCTION 

The objectives of our experiments were to measure the breakdown voltages at 
various pressures for the composite-structures of three different thicknesses and with 
three different diameter holes. (See Table 1.) The thickness refers to the thickness of the 
insulation between the copper layers for each composite structure on Figure 3. There 
were a total of twelve structures utilized in the experiments not including a few that were 
used to determine appropriate equipment set up. The twelve structures were organized 
into three groups of four. Each group consisted of a baseline structure without holes that 
was used to determine a baseline breakdown voltage for the operating pressures used. 
The other three consisted of one wafer each with 300, 400, 500pm diameter holes 
respectively. 


Group one 

127 pm structure 

Group two 

254 pm structure 

Group three 

381 pm structure 

Baseline 

300pm holes 
400pm holes 
500pm holes 

Baseline 

300pm holes 
400pm holes 
500pm holes 

Baseline 

300pm holes 
400pm holes 
500pm holes 


Table 1. Composite Structure’s Experimental Matrix. 


Each MSE was attached to leads from the high voltage DC power supply and 
placed inside the vacuum chamber which was subsequently evacuated to a minimum of 
10'^ Torr and then back filled with ultra pure research grade Argon (with a purity of 
99.995%) to the needed experimental pressures. After each measurement the vacuum 
chamber was once again evacuated to remove ionized Argon and any contaminants and 
then backfilled with Argon to the next pressure. The breakdown voltage was determined 
by increasing the applied voltage and monitoring the voltage with a Tektronix 
oscilloscope. Once breakdown voltage was achieved the voltage level was recorded. 
(See Appendix A for diagrams of equipment and operational procedures.) 
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Once the breakdown voltages were obtained for all microstruetures the data was 
put into MATLAB and graphed. (See appendix B for data tables of breakdown voltages 
and Pasehen eurves.) A typieal Paschen curve for Argon is given in Figure 3. The data 
was then analyzed to determine where the minimum breakdown voltage oeeurred for 
eaeh test case. 

B. CALCULATING EXPECTED RESULTS 

Once the experimental data was obtained, it was compared with the expected 
results. To determine the expected results a Pasehen curve was generated for eaeh 
strueture. The expected Paschen curves were found by determining the experimental 
eonstants using Equation (3): 


Cj 

ln(pJ) + ln(C2) 


(3) 


where C/ and C 2 are the experimental constants and Vb the breakdown voltage. Solving 
for the minimums using Equations (4) and (5) gives: 


P^min = 


2.718 

C 


_2.718*Ci 

' hmin 




(4) 

(5) 


where pdmin and Vb min are the minimum values from the experimental data. The 
following three tables (Tables 2, 3, 4) show the MATEAB caleulated C; and C 2 values 
for each micro structure: 


Table 2. 


Group 1 C, & C 2 Values 

Hole 

Diameter 

Ci (Volts’^) 

C 2 (Torr’kcm'^) 

Baseline 

SOOpm 

400|im 

500pm 

181820.00 

172910.00 

107800.00 

164040.00 

1764.90 

1926.20 

1153.60 

1741.60 

Group 1 C; & C 2 MAI 

XAB calculatec 


values. 
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Table 3. 


Table 4. 


Group 2 Ci &. C 2 Values 

Hole 

Diameter 

Ci (Volts’^) 

C 2 (Torr’kcm'^) 

Baseline 

71994.00 

627.1777 

300pm 

101420.00 

1111.6000 

400pm 

92142.00 

894.4320 

500pm 

95202.00 

1035.0000 

Group 2 Cl &C 2 MAI 

"LAB ealeulatec 


Group 3 Ci & C 2 Values 

Diameter 

C, (Volts'^) 

C 2 (Torr’^-cm'^) 

Baseline 

300pm 

400pm 

500pm 

50690.00 

60938.00 

92879.00 

38286.00 

488.5677 

707.8125 

934.9845 

369.0127 


values. 


Group 3 C] & C 2 MATLAB ealeulated 


values. 


Interesting to note is the C] and C 2 values deerease as the thiekness of the wafers 
inereases. The values also deerease as the hole diameter deereases. 

The next step was to graph the expeeted Pasehen eurve and eompare it to the 
eurve generated by the experimental data. 
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Paschen curve of thin wafer Baseline. 


600 


Paschen curve of thin wafer vwth 300um holes 




Paschen curve of thin wafer with 400um holes Paschen curve of thin wafer with SOOum holes 




Figure 6. Graphical Representation of group one (thin) wafer actual breakdown 
voltages and expected breakdown voltage curve. 
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Paschen cun« o1 middle wafer Baseline. Paschen curve of middle wafer with 300um holes 




Paschen cur\e of middle wafer with 400um holes Paschen curve of middle wafer with SOOum holes 



Figure 7. Graphical Representation of group two (middle) wafer actual breakdown 
voltages and expected breakdown voltage curve. 
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Paschen curve of thick wafer Baseline. Paschen cune of thick wafer with 300um holes 




Paschen cur\e of thick wafer with 400um holes Paschen curve of thick wafer with 500um holes 




Figure 8. Graphical Representation of group three (thick) wafer actual breakdown 
voltages and expected breakdown voltage curve. 


As can be seen in Figures 6, 7 and 8, the lowermost data points are close to the expected 
curves. Some variation can be seen but the minima are very close. 

C. REPEATABILITY 

The next question to be answered was whether the experiments can be duplicated 
with some degree of consistency. Two wafers were randomly chosen to be retested with 
the following results (Tables 5 and 6, Figures 9 and 10); 
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400iim - Group 2 structure with d=0.285mm/0.0112in 

First Run 

Second Run 

Pressure (p) (miiii-Torr) 

Breakdown Voitage (Volts) 

Pressure (p) (milli-Torr) 

Breakdown Voltage 
(Volts) 

11.2 

N/A 

N/A 

N/A 

30.3 

326 

30 

318 

51.3 

310 

50.7 

292 

83.7 

290 

80.9 

268 

102 

270 

101 

248 

150 

282 

150 

292 

203 

334 

200 

316 

498 

354 

499 

288 


Table 5. Data table eontrasting the first and seeond run with the same wafer. 


Paschen curve comparison of two runs of Group 3 wafer with 400um hoies 



400[im diameter holes. 


The eomparison between the first and seeond runs (as shown in figure 9) of the 

middle wafer with dOOpm diameter holes shows a high degree of similarity. The 

minimum voltages are within 22 volts of eaeh other at the same pressure. 
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500um - Group 3 structure with d=0.372mm/0.01465in 

First Run 

Second Run 

Pressure (p) (miiii-Torr) Breakdown Voitage (Voits) 

Pressure (p) (miiii-Torr) Breakdown Voitage (Voits) 

10.6 N/A 

30.7 346 

50.6 346 

83.1 308 

98.2 298 

149 298 

198 282 

498 352 

N/A N/A 

34 316 

50 306 

81.6 292 

N/A N/A 

150 296 

201 306 

N/A N/A 


Table 6. Data table contrasting the first and second run with the same wafer 


Paschen curve comparison of two runs of Group 3 wafer with 500um holes 



Figure 10. Comparison of Paschen curves for two runs of the thick wafer with 500pm 

diameter holes 


The comparison of the first and second runs (as shown in Figure 10) of the thick 
wafer with 500pm diameter holes shows that there is less consistency between these runs. 


22 
































The minimum voltage for breakdown varies by 10 volts but the pressure regions vary by 
116 Torr with the first run minima oeeurring at 198 Torr and the seeond run minima 
oeeurring at 81.6 Torr. 

The differenee between the two runs ean be in part eontributed to the wear on the 
eleetrodes of the mierostruetures, whieh modifies the three dimensional geometry 
slightly. 

D. DETERIORATION 

One of the main unknowns of this experiment is the amount of deterioration of the 
eleetrodes during the ionization proeess. The purpose of this experiment is not to 
quantify the amount of deterioration of the eleetrodes but it must be realized that the 
deterioration will ehange the three dimensionality of the geometry over time. The 
deterioration is a funetion of the time that the eleetrodes are under diseharge eonditions; 
the robustness of the material used for the eleetrodes; and the voltage-eurrent eonditions. 
In our ease eopper was used to demonstrate that deterioration did in faet take plaee on the 
MSE arrays that were examined both before and after the ionization proeess. Mieroseope 
images are presented below -please see Perry [7] for further details. 
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Figure 11. Photo of 300|^m diameter hole at 290 magnifieation before testing. 
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Figure 12. Photo of 400^m diameter hole at 48 magnification before testing. 
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Figure 13. Photo of SOO^m diameter hole at 48 magnification before testing 


In Figures 11,12 and 13 some roughness around the edges of the hole can be seen 
but these are due to the inaccuracies of the manufacturing methods used. 
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Figure 14. Photo of 300|.im diameter hole at 290 magnification after testing. 
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Figure 15. Photo of 400^m diameter hole at 48 magnification after testing. 
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Figure 16. Photo of SOO^m diameter hole at 290 magnification after testing 


As can be seen in Figures 14, 15 and 16, there is significant deterioration of the 
electrodes after the ionization process in Argon has occurred. The amount of 
deterioration varies depending on the hole size and the time that the hole was subjected to 
the discharge. In this respect, copper was not a robust electrode material - ion engines 
require long life times, therefore future research must include the study of suitable 
electrode materials. 
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V. CONCLUSIONS 


This thesis shows that the ionization of Argon, whieh is an alternate fuel to 
Xenon, ean be aehieved at low voltages by utilizing Miero-Struetured Eleetrode (MSE) 
arrays. The MSE arrays serve to eoneentrate the eleetrie fields between eleetrodes 
allowing the produetion of eleetrons enhaneing the field emission effeet. Minimum 
breakdown voltages between 240 and 280 volts at pressures of around 100 milli-Torr 
(0.133N/m ) were eonsistently obtained with arrays of hole diameter ranging from 300 to 
500pm. The MSE arrays were fabrieated from eomposite Cu-Si02-Cu struetures with 
varying thicknesses from 100 pm to 400 pm. The micro structures (holes) were 
fabricated using precision conventional machining and hole diameters between 300 and 
500 um were achieved. 

Using the experimental data and MATEAB we were able to obtain the 
experimental constants {Ci and C 2 ) for the Paschen curves. With these values in hand the 
minimum breakdown voltages were calculated for each of the MSE arrays. The values of 
the experimental constants decreased as the thickness of the layered structures increased 
and as the hole diameter decreased. 

Another achievement has been the utilization of commonly available materials 
and manufacturing methods for the construction of the MSE arrays. Thus yielding a 
low-cost method for conducting experiments in this field 
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VI. RECOMMENDATIONS 


Several reeommendations resulted from this study 

A. REPEATABILITY STUDIES 

Additional repeatability studies should be performed to verily that ionization ean 
be aehieved for the same MSE array many times at the same voltage. This will also serve 
to establish aeeurate empirieal data required for any future optimization studies. 

B. DETERIORATION STUDIES 

Further work should inelude the study of the rate and amount of deterioration that 
is ineurred by repeating the ionization proeess many times with the same MSE array. 
With eaeh repetition the three dimensionality of the struetures ehanged and the modeling 
and understanding of this ehange will be vital to predieting future performanee. 

C. FLOW STUDIES 

This experiment was eondueted in a statie atmosphere. Future ion engine 
applieations will not be utilizing sueh an atmosphere but will instead be utilizing the flow 
of Argon through the MSE array. Additional studies should be performed to determine 
the effeets of a flow (both eonstant and variable) through the array on the breakdown 
voltages and on ionization effieieney. 

D. MATERIAL STUDIES 

Another point of interest is the materials that make up the MSE array. For these 
experiments eopper eleetrodes were used. Ion propulsion applieations will require long 
life times, therefore suitable materials need to be identified 

E. STRUCTURES STUDIES 

The eurrent experiment used a eireular hole through the eleetrodes and insulator 
material. Other three dimensional struetures should be investigated for optimal ionization 
performanee and effieieney. 
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APPENDIX A 


A. PROCEDURE FOR VACUUM CHAMBER OPERATION. 

Steps: 

1. Vent vaeuum ehamber to normal atmospherie pressure by opening 
venting valve. 

2. Open vaeuum ebamber and eonneet wafer to eopper leads utilizing 
plastie elips. Ensure that tbe two eopper leads are attaebed on opposite 
sides of tbe wafer from eaeb other. 

3. Close vaeuum ehamber and inspeet seals. Ensure venting valve and 
roughing pump number 2 valve are elosed. Open isolation valve and 
start roughing pump number 1. 

4. Turn on instrument panel and thermoeouple (TC) gages 1 and 2 and 3. 
Onee TC2 and TC3 indieate ehamber pressure around lOOmilli-Torr turn 
power on to the Turbo Moleeular Pump (TMP) and monitor as power up 
sequenee eommenees. Turn on Eilament gage and monitor pressure. 

5. Onee vaeuum ehamber pressure is at its minimum pressure 
(approximately 1.5x10'^ Torr) fasten Argon line to vent valve on 
ehamber and ensure metering valve is elosed. Open pressure valve on 
Argon tank and open eut off valve to metering valve on Argon line to 
vaeuum ehamber. 

6. Close Isolation Valve and begin baekfilling vaeuum ehamber with Argon 
by opening vent valve and open metering valve. Adjust and eontrol flow 
of Argon into Vaeuum ehamber by adjusting metering valve and 
monitoring TC3 meter. 

7. Onee pressure in vaeuum ehamber has reaehed required pressure shut off 
metering valve and elose vent valve. 
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B. PROCEDURE FOR INSTRUMENT PANEL OPERATION 

Steps: 

1. Once the procedures in part A are complete and the experiment is ready to 
be run, turn on power to the equipment cabinet ensuring that the voltmeter, 
amp meter and oscilloscope power up. 

2. Switch amp meter from AC to DC setting. Volt meter should default to 
DC setting. Set oscilloscope to 5 volts per division and 5 seconds per 
division. 

3. For Safety Ensure No One is touching any wires or the vacuum chamber. 

4. Ensure power supply is set to zero volts and turn on power supply. 

5. Begin by incrementing voltage levels and monitor oscilloscope and 
voltmeter for breakdown voltage. 

6. Once breakdown voltage has been achieved record voltage and turn off 
power supply and set dials on power supply back to zero volts. 

7. Set up for next experiment. 
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DIAGRAMS 



Figure 17. Diagram of vacuum chamber assembly. 





Figure 18. Diagram of control panel for pressure sensors and pump controls. 
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^^^^ressme Valve 



Figure 20. Argon supply system layout. 





APPENDIX B 


A. DATA TABLES 
1. Group 1 


Baseline - Thin Wafer with d=0.154mm/0.00605in 

Pressure (p) (milli-Torr) 

Breakdown Voltage (Volts) 

13.5 

N/A 

32.4 

366 

51.4 

310 

82.3 

286 

100 

280 

153 

282 

200 

320 

495 

346 

Table 7. Group 1 baseline. 

300pm - Thin Wafer with d=0.137mm/0.00535in 

Pressure (p) (milli-Torr) 

Breakdown Voltage (Volts) 

13.8 

N/A 

35.2 

326 

51.3 

292 

82.2 

262 

103 

244 

151 

252 

202 

260 

509 

302 

Table 8. Group 1 with 300pm holes. 

400um - Thin Wafer with d=0.153mm/0.00605in 

Pressure (p) (milli-Torr) 

Breakdown Voltage (Volts) 

13.4 

N/A 

32.1 

364 

52.9 

310 

82.9 

276 

98.3 

264 

154 

254 

204 

276 

517 

348 


Table 9. Group 1 with 400pm holes. 


41 




500um - Thin Wafer with d=0.136mm/0.00535in 


Pressure (p) (miiii-Torr) Breakdown Voitage (Voits) 


13 

N/A 

35.1 

316 

50 

284 

83.8 

266 

102 

256 

148 

264 

199 

326 

500 

364 


Table 10. Group 1 with 500pm holes. 


2 . 


Group 2 


Baseiine - Middie Wafer with 
d=0.287mm/0.0113in 


Pressure (p) (miiii-Torr) Breakdown Voitage (Voits) 


12.9 

N/A 

31.5 

386 

50.4 

374 

85.9 

340 

100 

328 

151 

312 

198 

312 

496 

386 


Table 11. Group 2 baseline. 


300pm - Middie Wafer with d=0.303mm/0.0119in 


Pressure (p) (miiii-Torr)|Breakdown Voitage (Voits) 


11.9 

N/A 

33.9 

296 

50.9 

276 

80.7 

248 

106 

258 

159 

258 

208 

274 

503 

324 


Table 12. Group 2 with 300pm holes. 
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400pm - Middle Wafer with d= 

=0.284mm/0.01115in 

Pressure (p) (milli-Torr) Breakdown Voltage (Volts) 

10.4 

N/A 

31.8 

332 

48 

318 

84.3 

284 

107 

280 

152 

284 

200 

298 

500 

360 


Table 13. 

Group 2 with 400pm holes. 

500pm - Middle Wafer with d=0.26mm/0.01025in 

Pressure (p) (milli-Torr) Breakdown Voltage (Volts) 

11.9 

N/A 

31.7 

322 

56.6 

290 

81.4 

260 

101 

250 

150 

272 

200 

318 

499 

300 

800 

340 

Table 14. 

Group 2 with 500pm holes. 


3. Group 3 



Table 15. Group 3 baseline. 
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300pm - Thick Wafer with d=0.384mm/0.0151 in 

Pressure (p) (miiii-Torr) 

Breakdown Voitage (Voits) 

14.1 

N/A 

32.2 

382 

50.7 

320 

84.1 

284 

100 

234 

149 

262 

200 

276 

518 

348 

Table 16. Group 3 with 300pm holes. 

400pm - Thick Wafer with d=0.285mm/0.0112in 

Pressure (p) (miiii-Torr) 

Breakdown Voitage (Voits) 

11.2 

N/A 

30.3 

326 

51.3 

310 

83.7 

290 

102 

270 

150 

282 

203 

334 

498 

354 

Table 17. Group 3 with 400pm holes. 

500pm - Thick Wafer with d=0.372mm/0.01465in 

Pressure (p) (miiii-Torr) 

Breakdown Voitage (Voits) 

10.6 

N/A 

30.7 

346 

50.6 

346 

83.1 

308 

98.2 

298 

149 

298 

198 

282 

498 

352 


Table 18. Group 3 with 500pm holes. 
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B. PASCHEN CURVE RESULTS 

1. Paschen Curve Comparison of Data 


Paschen curves of thin (Group 1) wafers. 



Figure 21. Comparison of thin (Group 1) breakdown voltages. 
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Paschen curves of middle (Group 2} wafers. 



Figure 22. Comparison of middle (Group 2) breakdown voltages. 
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Paschen curves of thick (Group 3) wafers. 



Figure 23. Comparison of thick (Group 3) breakdown voltages. 
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APPENDIX C 


A. MATLAB CODE USED TO GRAPH FIGURES 6-8 

%% Thesis Data 

clear all 
close all 
clc 

%%%%Thinnest Wafers 

%% Thinnest wafer Baseline 

pdmin= 100*.0154 %milli-Torr 

pdmin=pdmin/1000 %Torr 

Vbmin=280 

c2=2.718/pdmin 

cl=Vbmin*c2/2.718 

x-00066:.0001:.01; 

vb=(c 1. *x) ./(log(x)+log(c2)); 

figure(l) 

subplot(2,2,l) 

plot(x,vb,'r.-Vlinewidth',3) 

hold on 

d=.0154; %cm or 0.00605in 

P=[32.4 51.4 82.3 100 153 200 495]; %% pressure in milliTorr 

P=P/1000; %% pressure in Torr 

Vb=[366 310 286 280 282 320 346]; %% voltage in volts 

pd=d*P; 

plot(pd,Vb,'linewidth',3) 

title('Paschen curve of thinnest wafer Baseline.') 

xlabelCpd (Torr-cm)') 

ylabel('Vb') 

legend('Expected','Data','Location','SouthEast') 
grid on 
hold on 
scatter(pd,Vb) 
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%% 300um thinnest wafer 

pdmin= 103*.0137 %milliTorr 

pdmin=pdmin/1000 %Torr 

Vbmin=244 

c2=2.718/pdmin 

cl=Vbmin*c2/2.718 

x-00066:.0001:.01; 

vb=(cl .*x)./(log(x)+log(c2)); 

subplot(2,2,2) 

plot(x,vb,'r.-Vlinewidth',3) 

hold on 

d=.0137; %cm or 0.00535in 

P=[35.2 51.3 82.2 103 151 202 509]; %% pressure in milliTorr 
P=P/1000; %% pressute in Torr 

Vb=[326 292 262 244 252 260 302]; %% voltage in volts 
pd=d*P; 

plot(pd,Vb,'linewidth',3) 

title('Paschen curve of thinnest wafer with 300um holes') 

xlabelCpd (Torr-cm)') 

ylabel('Vb') 

legend('Expected', 'Data', 'Location',' SouthEast') 

grid on 

%axis tight 

hold on 

scatter(pd,Vb) 

%% 400um thinnest wafer 

pdmin= 154*.0153 %milliTorr 

pdmin=pdmin/1000 %Torr 

Vbmin=254 

c2=2.718/pdmin 

cl=Vbmin*c2/2.718 

x-001:.0001:.01; 

vb=(cl .*x)./(log(x)+log(c2)); 

subplot(2,2,3) 

plot(x,vb,'r.-','linewidth',3) 

hold on 
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d=.0153; %cm or 0.00605in 

P=[32.1 52.9 82.9 98.3 154 204 517]; %% pressure in milliTorr 
P=P/1000; %% pressute in torr 

Vb=[364 310 276 264 254 276 348]; %% voltage in volts 
pd=d*P; 

plot(pd,Vb,'linewidth',3) 

title('Pasehen eurve of thinnest wafer with 400um holes') 

xlabel('pd (torr-em)') 

ylabel('Vb') 

legend('Expeeted', 'Data', 'Loeation',' SouthEast') 
grid on 
hold on 
seatter(pd,Vb) 

%% 500um thinnest wafer 

pdmin= 102*.0153 %millitorr 

pdmin=pdmin/1000 %torr 

Vbmin=256 

e2=2.718/pdmin 

el=Vbmin*e2/2.718 

x-00066:.0001:.01; 

vb=(e 1. *x) ./(log(x)+log(e2)); 

subplot(2,2,4) 

plot(x,vb,'r.-','linewidth',3) 

hold on 

d=.0153; %em or 0.00605in 

P=[35.1 50 83.8 102 148 199 500]; %% pressure in millitorr 
P=P/1000; %% pressute in torr 

Vb=[316 284 266 256 264 326 364]; %% voltage in volts 
pd=d*P; 

plot(pd,Vb,'linewidth',3) 

title('Pasehen eurve of thinnest wafer with 500um holes') 

xlabel('pd (torr-em)') 

ylabel('Vb') 

legend('Expeeted', 'Data', 'Loeation',' SouthEast') 
grid on 
hold on 
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scatter(pd,Vb) 


%%%%Middle Wafers 

%% Middle wafer Baseline 

pdmin= 151 *.0287 %millitorr 

pdmin=pdmin/1000 %torr 

Vbmin=312 

c2=2.718/pdmin 

cl=Vbmin*c2/2.718 

x-002:.0001:.015; 

vb=(c 1. *x) ./(log(x)+log(c2)); 

figure(2) 

subplot(2,2,l) 

plot(x,vb,'r.-Vlinewidth',3) 

hold on 

d-0287; %cm or 0.01 Din 

P=[31.5 50.4 85.9 100 151 198 496]; %% pressure in millitorr 
P=P/1000; %% pressute in torr 

Vb=[386 374 340 328 312 312 386]; %% voltage in volts 
pd=d*P; 

plot(pd,Vb,'linewidth',3) 

title('Paschen curve of middle wafer Baseline.') 

xlabel('pd (torr-cm)') 

ylabel('Vb') 

legend('Expected','Data','Location','NorthEast') 
grid on 
hold on 
scatter(pd,Vb) 

%% 300um middle wafer 

pdmin= 80.7*.0303 %millitorr 

pdmin=pdmin/1000 %torr 

Vbmin=248 

c2=2.718/pdmin 

cl=Vbmin*c2/2.718 

x-001:.0001:.016; 

vb=(c 1. *x) ./(log(x)+log(c2)); 
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subplot(2,2,2) 
plot(x,vb,'r.-Vlinewidth',3) 
hold on 

d-0303; %cm or0.0119in 

P=[33.9 50.9 80.7 106 159 208 503]; %% pressure in millitorr 
P=P/1000; %% pressute in torr 

Vb=[296 276 248 258 258 274 324]; %% voltage in volts 
pd=d*P; 

plot(pd,Vb,'linewidth',3) 

title('Pasehen eurve of middle wafer with 300um holes') 

xlabel('pd (torr-em)') 

ylabel('Vb') 

legend('Expeeted','Data','Loeation','NorthEast') 
grid on 
hold on 
seatter(pd,Vb) 

%% 400um middle wafer 

pdmin= 107*.0284 %millitorr 

pdmin=pdmin/1000 %torr 

Vbmin-280 

e2=2.718/pdmin 

el=Vbmin*e2/2.718 

x-003:.0001:.016; 

vb=(el .*x)./(log(x)+log(e2)); 

subplot(2,2,3) 

plot(x,vb,'r.-','linewidth',3) 

hold on 

d-0284;%em or0.01115in 

P=[31.8 48 84.3 107 152 200 500]; %% pressure in millitorr 
P=P/1000; %% pressute in torr 

Vb=[332 318 284 280 284 298 360]; %% voltage in volts 
pd=d*P; 

plot(pd,Vb,'linewidth',3) 

title('Pasehen eurve of middle wafer with 400um holes') 

xlabel('pd (torr-em)') 

ylabel('Vb') 
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legend('Expected', 'Data', 'Location',' SouthEast') 
grid on 
hold on 
scatter(pd,Vb) 

%% 500um middle wafer 

pdmin= 101 *.026 %millitorr 

pdmin=pdmin/1000 %torr 

Vbmin=250 

c2=2.718/pdmin 

cl=Vbmin*c2/2.718 

x-0015:.0001:.02; 

vb=(c 1. *x) ./(log(x)+log(c2)); 

subplot(2,2,4) 

plot(x,vb,'r.-','linewidth',3) 

hold on 

d=.026; %cm or 0.01025in 

P=[31.7 56.6 81.4 101 150 200 499 800]; %% pressure in milbtorr 
P=P/1000; %% pressute in torr 

Vb=[322 290 260 250 272 318 300 340]; %% voltage in volts 
pd=d*P; 

plot(pd,Vb,'linewidth',3) 

title('Paschen curve of middle wafer with 500um holes') 

xlabel('pd (torr-cm)') 

ylabel('Vb') 

legend('Expected', 'Data', 'Location',' SouthEast') 
grid on 
hold on 
scatter(pd,Vb) 

%%%%Thick Wafers 
%% Thick wafer Baseline 
pdmin= 152*.0366 %millitorr 
pdmin=pdmin/1000 %torr 
Vbmin=282 
c2=2.718/pdmin 
cl=Vbmin*c2/2.718 
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x-0025:.0001:.008; 

vb=(c 1. *x) ./(log(x)+log(c2)); 

figure(3) 

subplot(2,2,l) 

plot(x,vb,'r.-Vlinewidth',3) 

hold on 

d=.0366; %cm or 0.0144in 

P=[53.4 83.8 101 152 205]; %% pressure in millitorr 

P=P/1000; %% pressute in torr 

Vb=[338 314 304 282 324]; %% voltage in volts 

pd=d*P; 

plot(pd,Vb,'linewidth',3) 

title('Paschen curve of thick wafer Baseline.') 

xlabel('pd (torr-cm)') 

ylabel('Vb') 

legend('Expected','Data','Location','NorthEast') 
grid on 
hold on 
scatter(pd,Vb) 

%% 300um thick wafer 

pdmin= 100*.0384 %millitorr 

pdmin=pdmin/1000 %torr 

Vhmin-234 

c2=2.718/pdmin 

cl=Vbmin*c2/2.718 

x-0016:.0001:.025; 

vb=(c 1. *x) ./(log(x)+log(c2)); 

subplot(2,2,2) 

plot(x,vb,'r.-','linewidth',3) 

hold on 

d-0384; %cm or0.0151in 

P=[32.2 50.7 84.1 100 149 200 518]; %% pressure in millitorr 
P=P/1000; %% pressute in torr 

Vb=[382 320 284 234 262 276 348]; %% voltage in volts 
pd=d*P; 

plot(pd,Vb,'linewidth',3) 
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title('Paschen curve of thick wafer with 300um holes') 

xlabel('pd (torr-cm)') 

ylabel('Vb') 

legend('Expected', 'Data', 'Location',' SouthEast') 
grid on 
hold on 
scatter(pd,Vb) 

%% 400um thick wafer 

pdmin= 102*.0285 %millitorr 

pdmin=pdmin/1000 %torr 

Vhmin-270 

c2=2.718/pdmin 

cl=Vbmin*c2/2.718 

x-0015:.0001:.016; 

vb=(c 1. *x) ./(log(x)+log(c2)); 

subplot(2,2,3) 

plot(x,vb,'r.-','linewidth',3) 

hold on 

d=.0285; %cm or 0.0112in 

P=[30.3 51.3 83.7 102 150 203 498]; %% pressure in millitorr 
P=P/1000; %% pressute in torr 

Vb=[326 310 290 270 282 334 354]; %% voltage in volts 
pd=d*P; 

plot(pd,Vb,'linewidth',3) 

title('Paschen curve of thick wafer with 400um holes') 

xlabel('pd (torr-cm)') 

ylabel('Vb') 

legend('Expected', 'Data', 'Location',' SouthEast') 
grid on 
hold on 
scatter(pd,Vb) 

%% 500um thick wafer 
pdmin= 198*.0372 %millitorr 
pdmin=pdmin/1000 %torr 
Vbmin=282 
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c2=2.718/pdmin 

cl=Vbmiii*c2/2.718 

x=.004:.0001:.02; 

vb=(c 1. *x) ./(log(x)+log(c2)); 

subplot(2,2,4) 

plot(x,vb,'r.-Vlinewidth',3) 

hold on 

d-0372;%cm or0.01465in 

P=[30.7 50.6 83.1 98.2 149 198 498]; %% pressure in millitorr 
P=P/1000; %% pressute in torr 

Vb=[346 346 308 298 298 282 352]; %% voltage in volts 
pd=d*P; 

plot(pd,Vb,'linewidth',3) 

title('Paschen curve of thick wafer with 500um holes') 

xlabel('pd (torr-cm)') 

ylabel('Vb') 

legend('Expected', 'Data', 'Location',' SouthEast') 
grid on 
hold on 
scatter(pd,Vb) 

B. MATLAB CODE USED TO GRAPH FIGURES 21-23 IN APPENDIX B 

%% Thesis Data Comparison 

clear all 
close all 
clc 


%%%%Thinnest Wafers 
%% Thinnest wafer Baseline 
db=.0154; %cm or 0.00605in 

Pb=[32.4 51.4 82.3 100 153 200 495]; %% pressure in millitorr 

Pb=Pb/1000; %% pressute in torr 

Vbb=[366 310 286 280 282 320 346]; %% voltage in volts 

pdb=db*Pb; 
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%% 300um thinnest wafer 
d3=.0137; %cm or 0.00535in 

P3=[35.2 51.3 82.2 103 151 202 509]; %% pressure in millitorr 
P3=P3/1000; %% pressute in torr 

Vb3=[326 292 262 244 252 260 302]; %% voltage in volts 
pd3=d3*P3; 

%% 400um thinnest wafer 
d4=.0153; %cm or 0.00605in 

P4=[32.1 52.9 82.9 98.3 154 204 517]; %% pressure in millitorr 
P4=P4/1000; %% pressute in torr 

Vb4=[364 310 276 264 254 276 348]; %% voltage in volts 
pd4=d4*P4; 

%% 500um thinnest wafer 
d5=.0153; %cm or 0.00605in 

P5=[35.1 50 83.8 102 148 199 500]; %% pressure in millitorr 
P5=P5/1000; %% pressute in torr 

Vb5=[316 284 266 256 264 326 364]; %% voltage in volts 
pd5=d5*P5; 

%% Comparison plot for thinnest wafers 
figure(l) 

plot(pdb,Vbb,'b.-Vlinewidth',3) 
hold on 

plot(pd3, Vb3, 'go:'linewidth',2) 
hold on 

plot(pd4,Vb4,'rx-.','linewidth',2) 
hold on 

plot(pd5,Vb5,'k+-','linewidth',2) 
hold on 
grid on 

title('Paschen curves of thinnest wafers.') 

xlabel('pd (torr-cm)') 

ylabel('Vb') 

legend('Baseline', '3 OOum', '400um', '500um', 'Location', 'SouthEast') 

%%%%Middle Wafers 
%% Middle wafer Baseline 
db-0287;%cm or 0.01 Din 
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Pb=[31.5 50.4 85.9 100 151 198 496]; %% pressure in millitorr 

Pb=Pb/1000; %% pressute in torr 

Vbb=[386 374 340 328 312 312 386]; %% voltage in volts 

pdb=db*Pb; 

%% 300um Middle wafer 
d3-0303;%em or0.0119in 

P3=[33.9 50.9 80.7 106 159 208 503]; %% pressure in millitorr 
P3=P3/1000; %% pressute in torr 

Vb3=[296 276 248 258 258 274 324]; %% voltage in volts 
pd3=d3*P3; 

%% 400um Middle wafer 
d4=.0284; %em or 0.01115in 

P4=[31.8 48 84.3 107 152 200 500]; %% pressure in millitorr 
P4=P4/1000; %% pressute in torr 

Vb4=[332 318 284 280 284 298 360]; %% voltage in volts 
pd4=d4*P4; 

%% 500um Middle wafer 
d5=.026; %em or 0.01025in 

P5=[31.7 56.6 81.4 101 150 200 499 800]; %% pressure in millitorr 
P5=P5/1000; %% pressute in torr 

Vb5=[322 290 260 250 272 318 300 340]; %% voltage in volts 
pd5=d5*P5; 

%% Comparison plot for Middle wafers 
figure(2) 

plot(pdb,Vbb,'b.-Vlinewidth',3) 
hold on 

plot(pd3, Vb3, 'go:'linewidth',2) 
hold on 

plot(pd4,Vb4,'rx-.','linewidth',2) 
hold on 

plot(pd5,Vb5,'k+-','linewidth',2) 
hold on 
grid on 

title('Pasehen eurves of middle wafers.') 

xlabel('pd (torr-em)') 

ylabel('Vb') 

legend('Baseline', '3 OOum', '400um', '500um', 'Loeation', 'SouthEasf) 
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%%%%Thick Wafers 
%% Thick wafer Baseline 
db=.0366; %cm or 0.0144in 

Pb=[53.4 83.8 101 152 205]; %% pressure in milbtorr 
Pb=Pb/1000; %% pressute in torr 
Vbb=[338 314 304 282 324]; %% voltage in volts 
pdb=db*Pb; 

%% 300um Thick wafer 
d3-0384;%cm or0.0151in 

P3=[32.2 50.7 84.1 100 149 200 518]; %% pressure in millitorr 
P3=P3/1000; %% pressute in torr 

Vb3=[382 320 284 234 262 276 348]; %% voltage in volts 
pd3=d3*P3; 

%% 400um Thick wafer 
d4-0285;%cm or0.0112in 

P4=[30.3 51.3 83.7 102 150 203 498]; %% pressure in millitorr 
P4=P4/1000; %% pressute in torr 

Vb4=[326 310 290 270 282 334 354]; %% voltage in volts 
pd4=d4*P4; 

%% 500um Thick wafer 
d5=.0372; %cm or 0.01465in 

P5=[30.7 50.6 83.1 98.2 149 198 498]; %% pressure in millitorr 
P5=P5/1000; %% pressute in torr 

Vb5=[346 346 308 298 298 282 352]; %% voltage in volts 
pd5=d5*P5; 

%% Comparison plot for Thick wafers 
figure(3) 

plot(pdb,Vbb,'b.-Vlinewidth',3) 
hold on 

plot(pd3, Vb3, 'go:'linewidth',2) 
hold on 

plot(pd4,Vb4,'rx-.','linewidth',2) 
hold on 

plot(pd5,Vb5,'k+-','linewidth',2) 
hold on 
grid on 

title('Paschen curves of thick wafers.') 
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xlabel('pd (torr-cm)') 
ylabel('Vb') 

legend('Baseline', '3 OOum', '400um', '500um', 'Location'/SouthEast') 
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